A Physical Model for Baryons in Clusters of Galaxies 
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ABSTRACT 



> 

' The X-ray emission from clusters of galaxies is one of the best observational 

f^) . probe to investigate the distribution of dark matter at intermediate and high 

redshifts. Since the disposition of the intracluster plasma (ICP) responsible of 
the emission is crucial to link X-ray properties to the global properties of the 
dark matter halos, we propose a semi-analytical approach for the diffuse baryons. 
This comprises the following blocks: Monte Carlo "merging histories" to describe 
the dynamics of dark matter halos; the central hydrostatic disposition for the 
ICP; conditions of shock, or of closely adiabatic compression at the boundary 



with the external gas, preheated by stellar energy feedbacks. From our model we 
predict the L~T correlation, consistent with the data as for shape and scatter. 

1. Introduction 



^ ■ Groups and clusters of galaxies constitute cosmic structures sufficiently close to 

equilibrium and with sufficient density contrast (5^2 10 2 inside the virial radius R) 
as to yield definite observables. They are dominated by dark matter (hereafter DM), 
while the baryon fraction is observed to be less than 20%. The great majority of these 
baryons are in the form of diffuse plasma (ICP) with densities n ~ 10~ 3 cm -3 and 
virial temperatures kT ~ GM m^/10i? ~ 5 keV, and are responsible for powerful 
X-ray luminosities L ~ 10 44 erg/s by optically thin thermal bremsstrahlung. As the 
plasma is a good tracer of the potential wells, much better than member galaxies, the 
X-ray emission is a powerful tool to investigate the mass distribution out to moderate 
and high redshifts. The ICP temperature directly probes the height of the potential 
well, with the baryons in the role of mere tracers; on the other hand, the luminosity, 
with its strong dependence on density (L oc n 2 ), reliably probes the baryonic content 
and distribution. Statistically, a definite L-T correlation is observed (albeit with 
considerable scatter) , and this provides the crucial link to relate the X-ray luminosity 



functions with the statistics of the dark mass M or with that of the corresponding T. 

Many numerical experiments, using hydrodynamical N-body (see especially Gheller, 
Pantano & Moscardini 1997, Bryan & Norman 1997), provide a comprehensive tool 
to model the ICP distribution in the potential wells. However, such numerical ex- 
periments still do not have enough dynamic range to describe DM and ICP over the 
full range from ~ 50 Mpc associated with the large scale structures (which guide the 
ongoing mergers of DM halos), to the inner 50 kpc where the ICP yields a substantial 
contribution to L, and moreover do not include properly non-gravitational effects, 
which instead cannot be ignored. 

The simplest semi-analytical approach is constituted by the Self Similar model 
(Kaiser 1986) which include only gravity and assumes the ICP amount to be propor- 
tional to the DM at all z and M. This leads to a relation L oc T 2 conflicting with the 
observed correlation for rich clusters, which is close to L oc T 3 (David et al. 1993; 
Mushotzky & Scharf 1997). A missing ingredient is thought to be the stellar energy 
feedback by supernovae. 

The above features motivate us to develop a comprehensive semi-analytical model 
for the ICP which include the effects of stellar energy scales. First we assume that such 
energy input is efficient in depleting the potential wells of the clusters progenitors, 
at z — 1 -T- 2, and in pre- heating the intergalactic medium to temperatures in the 
range 7\ = 0.1 4- 0.8 keV as recently evidentiated in the outer cluster atmosphere 
(Henriksen & White 1996). Then we describe clusters evolution as a sequence of 
hierarchical merging episodes of the DM halos, associated in the ICP to shocks of 
various strengths (depending on the mass ratio of the merging clumps), which provide 
the boundary conditions for the ICP to re-adjust to a new hydrostatic equilibrium. In 
§3 we shall show that our predictions for the X-ray properties of clusters are consistent 
with recent data (see also Cavaliere, Menci & Tozzi 1997, hereafter CMT97). 

2. A physical model for the ICP 

The mass growth of dark halos leading to the formation of groups and clusters, 
develop through the accretion of diffuse matter or massive merging events between 
virialized halos. In our approach the history of such episodes is followed in the frame- 
work of the hierarchical clustering by Monte Carlo simulations. During the mass 
growth, the pre-heated ICP is recovered through shocks of variable intensity, de- 
pending on the temperature ratio between the accreted gas and the virialized plasma 
contained in the main progenitor. The ICP is reset to a new equilibrium after each 
episode of accretion or merging. 

In this framework it is possible to check the reliability of the assumption of equilib- 
rium. Roettiger, Stone & Mushotzky (1997) show that after a major merging event, 
i.e., with a mass ratio less than 2.5, the non thermal contribution to the pressure 
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Figure 1: Fraction of groups and clusters which have experienced at least one massive 
merging event with a mass ratio less then 2.5 in the last 2 Gyrs. Solid line: tilted 
CDM, dotted line: Open CDM, dashed line: ACDM (see N. Mend this volume). 

lasts less than two Gyrs. Adopting this as a conservative rule to identify disturbed 
objects, we can compute the fraction of clusters at redshift z = for which the hy- 
drodynamical equilibrium does not fully apply. As is shown in fig. [TJ this fraction 
is less than 20% at cluster scales in most FRW universes. Of course the situation is 
better for small objects, as they are in average older than rich clusters. 



2.1. Hydrodynamical equilibrium 

To compute the disposition of the ICP in equilibrium in a DM potential well, we 
have just to solve the hydrostatic equilibrium equation 
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n dr 
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where n(r) is the gas density profile, M(< r) is the mass contained within the radius 
r (dominated by DM) and P is the pressure. To solve the equation we need the 
boundary condition, i.e. the density of the gas at the virial radius n(R), and the 
state equation, that we write for a polytropic gas: 



P(r) 



kT(r) 
[irriH 



n 7 (r) 



(2) 



where m# is the proton mass, fi ~ 0.6 is the mean molecular weight of the ICP, and 
7, ranging from 1 to 5/3, is the polytropic index treated as a free parameter. The 
resulting profile is: 



n(r) = n{R) [l + (3^- — -) (<j>{R) - <f>(rj) 
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where (j) = V/o~ r is the adimensional gravitational potential and (3 = [imuOrjkT^^ 
with a r the line-of-sight velocity dispersion of the DM particles. Here T 2 is the 
temperature at the virial radius. The above equation can be considered a generalized 
/3-model (Cavaliere & Fusco-Femiano 1978), which reconduces to the usual isothermal 
case for 7 — > 1. Before computing the ICP distribution, we then need the boundary 
condition n(R). 



2.2. Physics of shocks 



The key boundary condition is provided by the dynamic stress balance P 2 = 
P\ + rriHnxvf, relating the exterior and interior pressures Pi and P\ to the inflow 
velocity v\ driven by the gravitational potential at the boundary. Here n\ is the 
baryon density external to the virial radius, and it is assumed to be unbiased respect 
to the universal value, i.e. n\ = QbP where p is the total matter density. We expect 
the inflowing gas to become supersonic close to R, when muv\> 2kT\. In fact, many 
hydrodynamical simulations of loose gas accretion into a cluster (from Perrenod 1980 
to Takizawa & Mineshige 1997) show shocks to form, to convert most of the bulk 
energy into thermal energy, and to expand slowly remaining close to the virial radius 
for some dynamical times. So we take R as the shock position, and focus on nearly 
static conditions inside, with the internal bulk velocity t> 2 << v i- 

The post-shock state is set by conservations across the shock not only of the 
stresses, but also of mass and momentum, as described by the Rankine-Hugoniot 
conditions (see Landau Sz Lifshitz 1959). These provide at the boundary the temper- 
ature jump T 2 /Ti, and the corresponding density jump g = n(R)/rii which reads 
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for a plasma with three degrees of freedom. Eq. f| includes both weak (with T 2 ~ T l5 
appropriate for small groups accreting preheated gas, or for rich clusters accreting 
comparable clumps), and strong shocks (appropriate to "cold inflow" as in rich clusters 
accreting small clumps and diffuse gas). From clusters to groups, the density jump 
g(T) lowers from the maximum value of 4 towards unity. 

Given the inflow velocity V\ and the nearly static post-shock condition t> 2 << 
it is possible to work out the explicit expression of the post-shock temperature T 2 in 
the form: 



kTo 



pLm H vlA\ + VT+~t) 2 7 3 e 2 



+ — e 



4 10 20 (1 + 



(5) 



where e = 15A;Ti/4/im//w 2 . For e ^> 1 the shock is weak and T 2 ~ T\ is recovered 
as expected. In the case of strong shocks, e € 1 and the approximation kT 2 ~ 
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Figure 2: a) The (3(T) parameter entering equation ^. b) Baryonic fraction respect 
to the universal value Qb for different values of 7. Continuous line: 7=1; dotted 
line: 7 = 1.1; dashed line: 7 = 1.2. 

—V(R)/3 + 3kTx/2 holds, where the second term is the contribution from the non- 
gravitational energy input. 

Now we can compute the ICP distribution, assuming a specific choiche for the 
potential well (in the following we use forms given by Navarro, Frenk & White 1996, 
but of course this is not mandatory). First, we compute /3(T), which decline from 
^ 1 for rich clusters, to ~ 0.4 for poor groups, where the stellar competes with the 
gravitational energy, as shown in fig. 0a. 

In fig. §a,b we show the emission-weighted temperature and density profile for two 
different values of 7 in the case of a rich cluster with M ~ 1O 15 M . A temperature 
profile with a mild decrease out to r ~ 1 Mpc is in agreement with the observations 
(Markevitch et al. 1997), pointing toward a value 7^ 1.2. The corresponding baryonic 
fraction lowers down with the mass scale by a factor of three from clusters to groups 
(see fig. 0b). 

3. The L-T correlation 

The X-ray luminosity of a cluster with temperature profile T(r) and density profile 
n(r) can be written: 

L^{9\T)) fdh^T^r). (6) 

In fact, before computing the L-T relation, the statistical effect of the merging 
histories has to be taken into account. For a cluster or group of a given mass (or 
temperature), the effective compression factor squared (g 2 ) is obtained upon aver- 
aging eq. [| over the sequence of the DM merging events; in such events, T 2 is the 
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Figure 3: Temperature and density profile for a cluster with M = 1O 15 M . Continu- 
ous line: 7 = 1; dashed line: 7 = 1.2. 

virial temperature of the receiving structure, and T\ is the higher between the stellar 
preheating temperature and that from "gravitational" heating, i.e., the virial value 
prevailing in the clumps being accreted. All that is accounted for in our model us- 
ing Monte Carlo simulations; these are based on merging trees corresponding to the 
excursion set approach of Bond et al. (1991), consistent with the Press & Schechter 
(1974) statistics (see CMT97). The averaged (g 2 ) is lower than the g 2 computed with 
a single temparature 7\, because in many events the accreted gas is at a tempera- 
ture higher than the preheating value. In addition, an intrinsic variance is generated, 
reflecting the variance intrinsic to the merging histories (see fig. |]a). 

The net result is shown in fig [|b. In agreement with the observations (David 
et al. 1993; Ponman et al. 1996), the shape of the average L — T relation flattens 
from L oc T 5 at the group scale (where the nuclear energy from stellar preheating 
competes with the gravitational energy) to L oc T 3 at the rich cluster scales. At larger 
temperatures the shape asymptotes to L oc T 2 , the self-similar scaling of pure gravity. 
Notice the intrinsic scatter due to the variance in the dynamical merging histories, 
but amplified by the n 2 dependence of L. We note that the shape of the L-T relation 
is little affected by change of 7. The average normalization formally rises like p 1 ^ 2 (z), 
where p is the effective external mass density which increases as (1 + z) 2 (Cavaliere & 
Menci 1997) in filamentary large scale structures hosting most groups and clusters. 
This implies, e.g., an increase of 30% at z = 0.3, consistent with the observations by 
Mushotzky & Scharf (1997). 

4. Conclusions 



The ICP state in the hierarchically evolving gravitational wells constitutes the 
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Figure 4: a) Tie factor (g 2 ) averaged over merging histories (solid line) with 96% con- 
fidence level (shaded region). The non-averaged g is plotted for comparison (dashed 
line), b) L-T relation compared with the data by David et al. (1993) for clusters 
(open squares) and Ponman et al. (1996) for groups (solid squares). 

focus of our new approach. We propose that such state follows suit, passing through 
a sequence of equilibrium condition that we compute semi-analytically. These com- 
putations comprise: the merging histories of the DM potential wells, obtained with 
a large statistics from Monte Carlo simulations of the hierarchical clustering; the 
inner hydrostatic equilibrium disposition, updated after each merging episode; and 
the boundary conditions provided by strong and weak shocks, or even by a closely 
adiabatic compression, depending on the ratio of the infall to the thermal energy in 
the preheated external medium. 

The results of our model depend on two parameters, the external temperature Ti 
and density n\, which are not free. Specifically, we use for T\ the range 0.1 4- 0.8 
keV provided by the literature on stellar preheating. The value of ri\ for rich clusters 
is related to the DM density by the universal baryonic fraction. The expression of 
the bolometric luminosity is proportional to g 2 = (^2,Mi) 2 , the square of the density 
jump at the bounding shock. The average of such factor over the merging histories 
coupled with (3{T) is what gives to the statistical L — T correlation the curved shape 
shown in fig. |]b. In addition, our approach predicts an intrinsic variance of dynamical 
origin due to the different merging histories, and built in the factor g 2 . Moreover the 
decreasing temperature profiles are in agreement with the published data and with 
the results from advanced simulations. A straightforward application is the prediction 
for the X-ray statistics in different FRW universes (see N. Menci, this volume, and 
Cavaliere, Menci & Tozzi 1998). 
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